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ABSTRACT 

Aims. We report on a 14 ks XMM-Newton observation of the hyperluminous infrared galaxy IRAS 09104+4109, which 
harbors a Fype 2 quasar in its nucleus. Our analysis was aimed at studying the properties of the absorbing matter and 
the Fe K complex at 6-7 keV in this source. 

Methods. We analyzed the spectroscopic data from the PN and the MOS cameras in the 0.4-10 keV band. We also 
used an archival BeppoSAX 1-50 keV observation of IRAS 09104+4109 to investigate possible variations of the quasar 
emission. 

Results. The X-ray emission in the EPIC band is dominated by the intra-cluster medium thermal emission. We found 
that the quasar contributes ~35% of the total flux in the 2-10 keV band. Both a transmission- (through a Compton- 
thin absorber with a Compton optical depth of re ~ 0.3, i.e. Nh ~5 x 10 23 cm -2 ) and a reflection-dominated (to > 
1) model provide an excellent fit to the quasar continuum emission. However, the value measured for the EW of Fe 
Ka emission line is only marginally consistent with the presence of a Compton-thick absorber in a reflection-dominated 
scenario, which had been suggested by a previous, marginal (i.e. 2.5a) detection with the hard X-ray (15-50 keV), 
non-imaging BeppoSAX/PDS instrument. Moreover, the value of luminosity in the 2-10 keV band measured by the 
transmission-dominated model is fully consistent with that expected on the basis of the bolometric luminosity of IRAS 
09104+4109. From the analysis of the XMM-Newton data we therefore suggest the possibility that the absorber along 
the line of sight to the nucleus of IRAS 09104+4109 is Compton-thin. Alternatively, the absorber column density could 
have changed from Compton-thick to -thin in the five years elapsed between the observations. If this is the case, then 
IRAS 09104+4109 is the first "changing-look" quasar ever detected. 

Key words. Galaxies: individual: IRAS 09104+4109 - Galaxies: active - Galaxies: nuclei - X-ray: galaxies 
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The X-ray observations of IRAS 09104+4109 collected igon with the BeppoSAX PDS flux, 101 concluded that the 

so far lend support to this suggestion. Franceschim et al. chandra spec trum is also reflection-dominated. 
(2000; F00) reported on the analysis of a BeppoSAX obser- 
vation of this source. They found that the X-ray spectrum 

below 10 keV is dominated by the intra-cluster medium 2. XMM-Newton observation and data reduction 

(ICM) thermal (fcT ~ 5.5 keV) emission. The detection of 

a weak signal in the 15-60 keV band with the PDS instru- IR-AS 09104+4109 was observed with XMM- 
ment was interpreted by F00 as the primary emission of the Newton (Jansen et al. 2001 and references therein) 
buried quasar, which emerges from an absorbing screen of on A P ril 27, 2003 for -14 ks. The EPIC PN and MOS ob- 
servations were carried out in the full frame mode using 
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Table 1. Best-fit spectral parameters of the EPIC spectrum. See Sect. [3] for details. Column: (1) spectral model; (2) 
temperature of the ICM "cool core" component (keV); (3) ICM metallicity; (4) 2-10 keV flux of the quasar component 
(10 -13 erg cm -2 s" 1 ); (5) 2-10 keV luminosity of the quasar component (10 44 erg s" 1 ); (6) column density of the absorber 
(10 23 cm" 2 ); (7) energy of the Fe Ka line (keV); (8) intensity of the Fe Ka line (10~ 6 ph/cm 2 /s); (9) EW of the Fe 
Ka line (eV); (10) reduced x 2 and number of degrees of freedom. 
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a (A) transmission model; (B) reflection model; b absorption-corrected line against absorption-corrected continuum; c with respect 
to the pure reflection component. 



tasks for processing the raw PN and MOS data files, 
respectively, to generate the relative linearized event files. 
X-ray events corresponding to patterns 0-12(0-4) for 
the MOS(PN) cameras were selected. Hot and bad pixels 
were removed. The event lists were furthermore filtered to 
ignore periods of high background flaring according to the 
method presented in Piconcelli et al. (2004) based on the 
cumulative distribution function of background lightcurve 
count-rates. Final net exposures of 10.8, 13.4, and 13.4 ks 
were obtained for PN, MOS1 and MOS2, respectively. The 
source photons were extracted for the PN(MOS) camera 
from a circular region with a radius of 37(40) arcsec, 
while the background counts were estimated from a larger 
(i.e. ~ 75 arcsec radius) source-free region on the same 
chip. Appropriate response and ancillary files for all the 
EPIC cameras were created using respectively RMFGEN 
and ARFGEN tasks in the SAS. Combined MOS1+MOS2 
spectrum and response matrix were created. 

Light-curves in the 0.2-2 keV and 2-15 keV band were 
extracted. Both light-curves are consistent with a constant 
flux level over the whole XMM-Newton exposure. 

3. Spectral modelling 

Both PN and MOS spectra were grouped to have a mini- 
mum of 20 counts per bin to allow the use of x 2 minimiza- 
tion technique and fitted simultaneously. Given the cur- 
rent calibration uncertainties, we discarded events below 
0.4 keV and above 10 keV. All fits were performed using 
the XSPEC package (vll.3) and included the Galactic col- 
umn density value of Nh=1-81 x 10 20 cm" 2 (Murphy et 
al. 1996). Best- fit parameter values are given in the source 
frame, unless otherwise specified. The quoted errors on the 
model parameters correspond to a 90% confidence level for 
one interesting parameter (A% 2 = 2.71; Avni 1976). A cos- 
mology with (AmjOa) = (0.3,0.7) and a H — 70 km s" 1 
Mpc -1 is assumed throughout. 

The high-resolution ( £ 1 arcsec) Chandra observa- 
tion of IRAS 09104+4109 presented by 101 has defini- 
tively established that the emission from the ICM domi- 
nates the X-ray spectrum below 10 keV. Furthermore IRAS 
09104+4109 is a massive cooling-flow cluster (Fabian & 
Crawford 1995) showing a large radial temperature gradi- 
ent, from fcT w 3.3 keV in the bright, cool core up to fcT w 
7.8 keV at a distance of 200 kpc, i.e. 36 arcsec (101). For this 
reason, we initially fitted the combined PN+ MOS spectrum 
with a model consisting of two thermal plasma emission 



components (i.e. MEKAL in XSPEC) plus an absorbed power 
law. We fixed the temperature of one thermal component 
to 7.8 keV in order not to underestimate the total ICM 
emission in the hard X-ray band. Since the photon index of 
the power law was loosely constrained (r f»1.4±l) we fixed 
it to 1.8, which is the average value typically observed for 
radio-quiet quasars (e.g. Piconcelli et al. 2005). 

This model (transmission model hereafter) gave 
an excellent description of the EPIC data with a 
X 2 (dof) =0.94(219) (sec Fig.HJ. The best-fit values of the 
spectral parameters arc listed in Table 1. The value of ~ 
4.8 x 10 23 cm" 2 found for the column density translates in 
a Compton optical depth tc ~ 0.3 of the absorbing screen. 
The value of temperature (fcT = 3.9 j^'g keV) in the cluster 

core and metallicity (Z/Z 3 = 0A7±%\\%) of the ICM are 
consistent with 101. 

Works based on BeppoSAX and Chandra observa- 
tions favored an interpretation of the spectrum of IRAS 
09104+4109 below 10 keV in terms of reflection-dominated 
emission. We therefore replaced the absorbed power law 
in the transmission model with a Compton reflection com- 
ponent from neutral matter (i.e. PEXRAV model in XSPEC). 
For this spectral component, which is due to the repro- 
cessing of the emission from the obscured primary X-ray 
source, we assumed a T = 1.8 for the photon index of the 
incident power law, along with an inclination angle of i = 
50 deg (Tran et al. 2000) and solar metallicity for the re- 
flector. This model (reflection model hereafter) yielded an 
equally good fit to the XMM-Newton spectrum with a final 
X 2 (dof)= 0.96(220) (see Table 1 and Fig.©. According to 
both spectral models, the contribution of the quasar com- 
ponent to the total flux in the 2-10 keV band is «30-35%. 

As shown in Figs. [T^ and there is a prominent line- 
like emission feature at 4-5 keV (observer-frame) which is 
broader than the instrumental resolution at this energy and 
most likely due to a blend of lines associated with Fe K 
emission. This complex is partly accounted for by the strong 
FeXXV Ka emission line at 6.7 keV and the FeXXVI Ka 
emission line at 6.97 keV (which should be likely blended 
with a weak Fe K/3 line at 7.06 keV due to the reprocessing 
of the quasar continuum) from the two-temperature ICM. 
However, significant (at >99% confidence level) positive 
residuals are still present. We modelled this excess with an 
unresolved Gaussian emission line at 6.38±g;g| keV. This 
energy is consistent with a range of ionization states from 
Fel to FeXVI (Kallman et al. 2004), as typically observed 
in quasars (Jimenez-Bailon et al. 2005). In the "reflection- 
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Fig. 1. (a) Left: XMM-Newton PN (top) and MOS (bottom) spectra of IRAS 09104+4109 when the transmission model 
is applied. The lower panel shows the deviations of the observed data from the model in unit of standard deviations, (b) 
Right: Best-fit model for the transmission scenario. The different spectral components are also plotted (e.g. Section [3] for 
further details). 




dominated" scenario we measure an equivalent width of the 
Fe Ka line at ~6.4 keV of EWjf a = 402±^ eV (calculated 
with respect to the Compton reflection component). For the 
transmission model (i.e. assuming that the line and the con- 
tinuum are both absorbed) we derive EWk a = 390+212 e V 
(see Table 1). 



3.1. A comparison with broadband 1-50 keV BeppoSAX data 

To investigate on the possible year-timescale variability of 
the overall continuum spectral shape, we plot in Fig. [3^ 
the unfolded 2003 PN and 1998 BeppoSAX MECS+PDS 
spectra. These data have been unfolded through the in- 
strument response with respect to the best-fit model found 
by F00 (i.e. ICM thermal emission component + absorbed 
powerlaw + reflection + narrow Gaussian line at 6.4 keV; 
model sax98 hereafter). As expected, the 1-10 keV XMM- 
Newton and BeppoSAX spectra have similar shape being 
dominated by the ICM emission. The different normaliza- 
tions of the two spectra can be ascribed to the different 
source extraction regions. In particular, a fraction of ICM 
emission from the outskirt region of the cluster spreads out- 
side the PN/ MOS extraction radius. Interestingly, the PDS 
data above 10 keV appear to lie slightly above the extrap- 
olation of the MECS/ PN continuum level. We calculated 
a 20-30 keV flux level of F 2 o_ 3 o ~ 6.1 x 10~ 13 and ~ 



1.3 x 10 12 erg cm 2 s 1 for the transmission and re- 
flection model, respectively. The 20-30 keV flux of IRAS 
09104+4109 measured by BeppoSAX is F 20 - 3 o = 2.55±i;^ 
x 10~ 12 erg cm~ 2 s _1 , whereby the extrapolated XMM- 
Newton flux in this band is fainter by a factor of ~1.6-7.3, 
if the transmission model is assumed. On the other hand, 
the 20-30 keV flux estimated by reflection model is consis- 
tent with the BeppoSAX value. 

Fig.[3jD shows the iso-% 2 contour plot of the intensity of 
the Fe Ka emission line at 6.4 keV (Iko) versus the inten- 
sity of the reflection continuum for the 1998 BeppoSAX and 
2003 XMM-Newton observations. The XMM-Newton val- 
ues were derived by the reflection model, while for the 
BeppoSAX data we employed the sax98 model (tempera- 
ture and abundance of the ICM were permitted to vary 
within the 90% confidence interval measured for these pa- 
rameters with XMM-Newton, e.g. Table 1). These measure- 
ments taken at different epochs are consistent within 3rj 
errors, as expected if Compton hump and cold iron emis- 
sion arise from distant material. Furthermore, the XMM- 
Newton measurements are also compatible with the values 
of l Ka = 7.2±f l x 10~ 6 photons cm" 2 s" 1 and EW Ka 

1 In the case of this model, we also added to the model a 
Compton reflection component with R = 1 (where R is the solid 
angle in units of 2n subtended by the reflecting material) as 
commonly observed in the 10-50 keV AGN spectra (e.g. Risaliti 
2002; Reeves et al. 2006). 
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Fig. 3. (a) Left: Unfolded spectra for the 2003 XMM-Newton and 1998 BeppoSAX observations, (b) Right: Confidence 
contour plot showing the intensity of the Fe Ka emission line (in units of ph/cm 2 /s) against the intensity (i.e. pho- 
ton flux at 1 keV) of the reflection continuum (in units of ph/keV/cm 2 /s) obtained using the XMM-Newton and 
BeppoSAX (MECS+PDS) data. The contours are at 68%, 90% and 99% confidence levels for two interesting parameters. 



= l.lloia keV reported by 101 on the basis of a 1999 
Chandra observation of IRAS 09104+4109. The mean val- 
ues of EWifc were significantly larger, but the errors were 
also very large in the low S/N BeppoSAX and Chandra spec- 
tra. 



4. Discussion 

The XMM-Newton observation presented here has con- 
firmed the presence of a heavily absorbed quasar in 
the nucleus of the hyperluminous infrared galaxy IRAS 
09104+4109. The EPIC energy range is dominated by the 
thermal emission of the ICM and we estimated that the 
quasar contributes approximately ~30-35% of the total 2— 
10 keV flux. We were able to accurately study the Fe K 
complex at —6-7 keV in this source, modelling it with three 
distinct narrow emission lines from Fel-XVI Ka, FeXXV 
Ka and FeXXVI Ka transition, respectively The highly 
ionized Fe K lines originate in the ICM, while the line from 
cold iron is likely due to reprocessing of the quasar contin- 
uum off the circumnuclear environment. The better quality 
of EPIC data has allowed to overcome the problems of lim- 
ited spectral resolution and statistics of BeppoSAX data 
(e.g. F00), which prevented a correct estimate of the con- 
tribution from the different ions to the Fe emission complex. 

The quasar emission can be equally well fitted by either 
a transmission or a reflection model. The former implies a 
Compton-thin absorber with Nh ~ 5 x 10 23 cm" 2 , while 
the latter suggests a scenario where the primary X-ray 
continuum is blocked by a Compton-thick (tq ^ 1, i.e. 
Nh — c^ 1 ~l-6 x 10 24 cm" 2 ) obscuring screen and the 
emission observed by XMM-Newton is due to indirect 
radiation scattered into the line of sight (it is generally 
assumed that the absorber and the reflector are the same 
material, i.e. the torus, e.g. Matt et al. 1996; Molendi et al. 
2003). Using the present data these models are statistically 
indistinguishable. However, an EW^a ~ 400±200 eV mea- 
sured assuming the reflection model is significantly lower 
than expected for a truly "reflection-dominated" spectrum 
with [Fe/H] = 0. In fact, the most prominent feature in the 
2-10 keV spectrum of a heavily obscured AGN, as in the 
case of IRAS 09104+4109 (2 x 10 24 < N H £ 10 25 cm" 2 , 
e.g. 101), is a strong Fe Ka emission line with an EWku 



1 keV (Levenson et al. 2002; Guainazzi et al. 2005). 
This solid observational evidence matches well with the 
EW value predicted by the theoretical calculations (e.g. 
Leahy & Creighton 1993; Ghisellini et al. 1994; Matt et 
al. 1996). However, given the uncertainty in the reflection 
continuum flux due to the presence of the diffuse cluster 
emission, we also conservatively estimated the EW^q by 
using the lowest (highest) value in the 90% confidence 
interval measured for the normalization of the reflection 
component (Fe Ka emission line) and vice versa. We 
obtained a range of EW_r- q values spanning from —165 to 
870 eV, which is marginally consistent with the EWko ~ 
1 keV expected for a reflection-dominated scenario. As the 
continuum flux estimated for the obscured quasar with the 
reflection model (F 2 _i = 4.68±^ x 10" 13 erg cm" 2 s" 1 ) 
is similar to the value measured with the spatially resolved 
Chandra spectrum (F2-10 ~ 3.9 x 10~ 13 erg cm" 2 s" 1 , 
e.g. 101), the marginal discrepancy in the EWjf Q values 
between the two observations arises from the different 
best-fit line flux intensity. 

On the other hand, the EW^ Q = 390^^2 eV ( note the 
large error bars) found in the case of the transmission model 
is in agreement with the expected EW value of a Fe Ka line 
transmitted through an absorbing screen with Nh ~5 x 
10 23 cm" 2 is (-200-400 eV; e.g. Awaki et al. 1991; Leahy 
& Creighton 1993; Ghisellini et al. 1994). A scenario with 
a Compton-thin absorber along the line of sight to the nu- 
cleus of IRAS 09104+4109 is therefore physically plausible. 
Nonetheless, the flux difference by a factor of —1.6-7.3 be- 
tween the 2003 XMM-Newton and 1998 BeppoSAX obser- 
vation (a time interval of 3.5 yr at the source frame after 
the correction for the time dilation due to the z = 0.442) in 
the 20-30 keV band is puzzling. In fact, given the observed 
power spectral densities of low black hole mass (i.e. Mbh ~ 
10 7 Mq) Seyfert galaxies (Markowitz et al. 2003) and scal- 
ing linearly with the black hole mass (assuming Mbh = 2.4 
x 10 9 M for the nuclear black hole in IRAS 09104+41093), 
flux variations of a factor of ^ 2 in 3.5 yr should be con- 



We estimated a mass of 2.4 x 10 9 M for the black hole 
in the nucleus of IRAS 09104+4109 using the measurement of 
the Mgll line width from Hines & Wills (1993), i.e. FWHM = 
10,000 km/s, and the formula in Willott et al. (2003). 
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sidercd unlikely (e.g. Fiore et al. 1998). However, the major 
problem in the interpretation of this mismatch concerns the 
accuracy of BeppoSAX/PDS measurement. According to 
our analysis, the signal-to-noise ratio of these data is poor 
(only — 2.5c between 15 and 50 keV). This implies that any 
conclusion based only on the PDS data about the Compton- 
thick or Compton-thin nature of the absorber in the nucleus 
of IRAS 09104+4109 could be misleading. The addition of 
a primary continuum power-law component modified by a 
Compton-thick absorber (i.e. similarly to the best-fit model 
proposed by F00 for the BeppoSAX broadband spectrum) 
might easily account for the above flux discrepancy, but, as 
stressed above, a scenario with a Compton-thick absorber 
for the EPIC spectrum is at odds with the inferred quite 
small value (i.e. < 1 keV) of the Fe Ka EW. 

Given the large 1.3° FWHM PDS field of view (FOV), 
PDS data might be affected by the contamination of 
very hard X-ray sources possibly located outside the 
XMM-Newton and Chandra 15 arcmin radius FOVs. If a 
power- law source with T = 1.8 is assumed, the 20-30 keV 
flux measured for the PDS source (-F20-30 —2.55 x 10~ 12 
erg cm~ 2 s _1 ) is translated to -F^-io^ 7.5 x 10~ 12 erg 
cm~ 2 s _1 . According to the RXTE all-sky slew survey log 
N - log S function in the 2-10 keV band (Revnivtsev et 
al. 2004), at this flux level 0.025 sources are expected in a 
1.5 deg 2 area, which translates into a probability of 2.4%. 
Moreover, using the NED and SIMBAD catalogs we found 
a probable contaminating source in the H2O maser galaxy 
NGC 2782 (z = 0.008), which is located at -50 arcmin 
away from IRAS 09104+4109. The Chandra spectrum 
of this source (with a nuclear X-ray flux of a few times 
10~ 13 erg cm~ 2 s _1 ) strongly suggests the presence of 
a heavily (N H > 10 24 cm" 2 ) absorbed AGN (Zhang et 
al. 2006), whereby it is likely that NGC 2782 provides 
a sizable contribution to the 20-30 keV flux measured 
with the PDS. Assuming a pure reflection model, we 
estimated a F20-30 ^6 x 10~ 13 erg cm~ 2 s _1 , which must 
be considered a lower limit of the 20-30 keV emission 
from NGC 2782, because of the likely presence of the 
nuclear continuum emerging after transmission through 
the absorber. 

Alternatively, IRAS 09104+4109 could represent the 
first example of a "changing-look" quasar (Guainazzi et 
al. 2002; Matt et al. 2003 and reference therein) ever de- 
tected given that the XMM-Newton data are better ex- 
plained by a transmission-dominated model with Nh —5 
x 10 23 cm~ 2 while, if the PDS emission is entirely due 
to IRAS 09104+4109, the BeppoSAX data are consistent 
with a reflection-dominated spectral state. In fact, similar 
spectral transitions from a Compton-thick to a Compton- 
thin state (or vice-versa) have been observed, but only in 
Seyfert-like AGNs so far. This scenario implies a dramatic 
decrease (a factor of £ 5-10) in the line-of-sight absorbing 
column density during a timescale of 5 years and, in turn, 
suggests the presence of a largely inhomogeneous obscur- 
ing circumnuclear gas (Elvis et al. 2004; Elitzur & Shlosman 
2006). We estimated the line-of-sight crossing-time of an ob- 
scuring cloud in Keplerian motion around the central black 
hole of 2.4 x 10 9 M Q to explain the possible transition 
from a reflection-dominated to a transmission-dominated 
spectrum. We assumed a scenario similar to that described 
for NGC 1365 in Risaliti et al. (2007) (or NGC 3227, e.g. 
Lamer et al. 2003) where they found that the Compton- 



thick obscuring material responsible of the spectral tran- 
sition is located in the broad line region (BLR) and the 
size of the X-ray emitting region is £ 100 Rg. We esti- 
mated a distance of the BLR Rblr — 0.14 pc using the 
empirical relation Rblr-Mbh in Kaspi et al. (2000). The 
redshift-corrected crossing-time (Guainazzi et al. 2002) of 
a Keplerian cloud covering a region of size 100 Rg around 
a black hole of 2.4 x 10 9 M Q is —1.7 yr. Such a value is 
therefore consistent with the 3.5 yr (source frame) elapsed 
between the BeppoSAX and XMM-Newton observation. 

Finally, we calculated a ratio rx,boi = ^2-iofcev/Lboi 
= 0.016 using the 2-10 keV luminosity of L2-wkeV = 7.95 
x 10 erg s _1 measured for the transmission-scenario 
(see Table 1), and a bolometric luminosity of Lboi(— Lir) 
= 4.7 x 10 46 erg s _1 , which is largely dominated by the 
obscured quasar (e.g. Hines et al. 1999; Spoon et al. 2007). 
We also derived the value of rxfiol expected for IRAS 
09104+4109 using the value of W„(2500A)/L bo i typical for 
quasars reported in Elvis et al. (1994), correcting the Lboi 
value by a factor of 1 /3 as suggested by Fabian & Iwasawa 
(1999) in order not to count twice the UV emission, and 
the spectral index aox between 2500A and 2 keV, defined 
as aox = -0.384^(2 keV)/Z ! ,(2500A)] (e.g. Tananbaum 
et al. 1979). In particular, we used the relation aox 
= 0.137%(Z„(2500A)) - 2.638 reported by Steffen et al 
(2006). We converted from the monochromatic value of 
vl v (2 keV) to the L 2 -iokeV value by multiplying by a factor 
of 1.61, applying a photon index of T = 2. We obtained 
that r x ,boi = L 2 -wkeV /Lboi = 0.043 x L bo i,45 _0 ' 357 (with 
Lboi,45=Lboi /10 45 erg s _1 ). This implies that the expected 
value of r x ,boi for IRAS 09104+4109 is 0.011, which is 
close to the value of rx.boi — 0.016 measured using the 
derived for the transmission scenario. Since the 
observed luminosity of L,2-wkeV = 2.05 x 10 44 erg s _1 in 
the reflection model should be just few percent of the 
intrinsic one (Comastri 2004; 101), this result lends further 
support to the hypothesis of a Compton-thin absorber 
along the line of sight to the nucleus of IRAS 09104+4109. 



5. Summary 

The XMM-Newton observation of IRAS 09104+4109 sug- 
gests the possibility that the absorber along the line of 
sight to the nucleus of IRAS 09104+4109 is Compton-thin. 
If this is the case, it implies a scenario completely dif- 
ferent from that reported for this Type 2 quasar so far. 
It is worth stressing, however, that previous X-ray stud- 
ies of IRAS 09104+4109 inferred a reflection-dominated 
nature of its 2-10 keV spectrum mainly on the basis of 
poor signal-to-noise ~ 15-50 keV data taken with the non- 
imaging BeppoSAX/PDS detector. Future imaging spec- 
troscopy of IRAS 09104+4109 performed in the 10-50 keV 
range, say with Simbol-X or XEUS, is therefore needed to 
make definitive progress in measuring the exact continuum 
emission from the quasar and constraining the column den- 
sity of the nuclear absorber. A deep Chandra observation of 
IRAS 09104+4109 would also be useful to accurately quan- 
tify the strength of the Fe Ka emission line at 6.4 keV, 
which is a proxy for the Compton thickness of the absorber. 
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